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We consider mesons composed of light and heavy quarks and discuss the construction of the cor- 
responding meson wave functions in soft- wall AdS/QCD. We specifically take care that constraints 
imposed by chiral symmetry breaking and by the heavy quark limit are fulfilled. The main results 
are i) the wave functions of light mesons have a nontrivial dependence on the current quark mass, 
which gives rise to a mass spectrum consistent with the one including explicit breaking of chiral 
symmetry; ii) the wave functions of heavy-light mesons generate their correct mass spectrum, the 
mass splittings of vector and pseudoscalar states, and the correct scaling of leptonic decay constants 
fgq ~ I/^/ttTq; iii) the wave functions of heavy quarkonia produce their correct mass spectrum 
and lead to a scaling behavior of the leptonic decay constants /qq ~ y/ m Q and / c j ~ m c /V m & at 
m c <C m,b, consistent with potential models and QCD sum rules. 
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The last decade has been marked by significant prog- 
ress in the development of AdS/QCD - a new class of 
approaches based on gauge/gravity duality which try 
to model strong interactions in terms of fields propagat- 
ing in extradimensional curved manifolds, e.g. in anti-de 
Sitter (AdS) space. Fields in AdS^+i space are classified 
by unitary, irreducible representations of the SO(d, 2) 
group @, El| , which is isomorphic to the conformal group 
acting on the boundary of AdS space, where the dual 
conformal field theory (CFT) is living. With some as- 
sumptions the CFT can be truncated to QCD, and there- 
fore the AdS fields could be holographically matched to 
QCD operators and bound states. There are two main 
types of AdS/QCD approaches: top-down (brane con- 
structions in string theories leading to low-energy gauge 
theories with properties of QCD) and bottom-up models 
(phenomenological frameworks specifying the geometry 
of AdS space and bulk fields in order to incorporate the 
basic properties of QCD). In this paper we focus on one 
of the successful examples of bottom-up approaches — 
the soft- wall model (see e.g. Refs. @-[6|), which is based 
on a soft breaking of conformal invariance via the intro- 
duction of a dilaton field (holographic analogue of the 
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gluon condensate Q) in the exponential prefactor or the 
effective potential. It leads to a truncation of AdS space 
in the infrared and therefore provides confinement of bulk 
fields, which are expanded in a tower of massive Kaluza- 
Klein modes identified with radial excitations of hadrons. 
The dilaton field is also responsible for the mechanism of 
spontaneous breaking of chiral symmetry. 

The profiles of bulk fields in extra dimension are 
matched to hadronic wave functions. It has been shown 
that the extradimensional coordinate can be identified 
with the transverse impact variable characterizing the 
separation of partons in a hadron in light-front QCD 
(Light-Front Holography) Q. It is essential to use the 
quadratic profile for the dilaton field to generate Regge 
trajectories for hadron masses and to reproduce the cor- 
rect scaling of hadronic form factors at large values of 
Euclidean transverse momentum squared. It is also help- 
ful to use this profile since most of the calculations can 
be performed analytically. Of course, not all features 
of QCD have been yet incorporated in the formalism of 
the soft-wall model. The work on providing an accurate 
matching to QCD is in progress. Therefore, some objec- 
tions concerning the soft-wall model seem premature. 

This work is addressed to the problem of constructing 
the hadronic wave functions using the AdS and light- 
front QCD correspondence. Originally the idea of such 
correspondence was proposed in Refs. [9(. It was shown 
that, from the matching of matrix elements for physical 
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processes (e.g. from the electromagnetic or gravitational 
form factors of hadrons), one can relate the string mode 
- the bulk profile of the AdS field in a holographic di- 
mension, and the transverse part of the hadronic light- 
front wave function (LFWF) for the case of massless 
quarks. Later in Ref. [Io| |. in the case of a two-parton 
state the LFWF was generalized by the explicit inclu- 
sion of the constituent quark masses in the LF kinetic 
energy Y^iO^Xi + m i)/ x i- I R the LFWF this corresponds 
to the introduction of the longitudinal WF corresponding 
to the so-called Brodsky-Huang-Lepage (BHL) or Gaus- 
sian ansa tz flU ] (see also discussion in Ref. [Hj]). In 
Refs. [3, Il4| we studied the problem of the longitudi- 
nal part of the LFWF. In particular, in 13| following 
the ideas of Ref. , we derived the longitudinal part 
of the LFWF using constraints of heavy quark effective 
theory. It was based on the BHL ansatz, where the di- 
mensional parameter depends on the flavor of the con- 
stituent quarks. Note that the BHL ansatz deals with 
constituent quarks while the direct parameters of QCD 
are the current quarks. The main objective of this pa- 
per is to construct the longitudinal part of the LFWF 
in terms of current quark masses instead of constituent 
ones. In particular, we will show that this construction 
helps to introduce the mechanism of explicit breaking of 
chiral symmetry. This means that the explicit breaking 
of chiral symmetry is a property of the longitudinal part 
of the hadronic LFWF. The idea that explicit breaking 
of chiral symmetry is induced via the current quark mass 
dependence of the longitudinal LWFW was proposed in 
two-dimensional large N c QCD [l5| . This mechanism was 
later used in the context of the two-dimensional massive 
Schwinger model I16l- 18II. Recently this problem was re- 
examined in Refs. [19|, |20j. 

In our approach the breaking of the conformal and chi- 
ral symmetries is related to the presence of the dilaton 
field. We show that both symmetries are broken dynam- 
ically due to the interaction with a dilaton - background 
field living in the z direction of AdS space; in other words, 
it does not show up at the AdS boundary. The dilaton is 
a massless zero mode. As it was mentioned before we use 
the quadratic dilaton profile, which can be considered as 
the expectation value of the scalar bulk field with dimen- 
sion 2, which is holographically dual to the dimension-2 
gluon operator Ojp = (^, m in) plj - In different types of 
AdS/QCD approaches a few scenarios of chiral symmetry 
breaking have been proposed due to the presence of spe- 
cific background or scalar fields. These fields can be con- 
sidered as duals to the dimension-4 gluon Oa? = tr(G 2 ) 
and the dimension-3 quark Oq q = qq operators. Their 
couplings to the bulk fields integrated over z can be in- 
terpreted as holographic analogues of in-hadron conden- 
sates of the operators Oa^, Og 2 an d 0«j 9 - The reason 
for this assignment is the following: in AdS/QCD (espe- 
cially in the soft-wall model) matrix elements are defined 
by integrals over profiles of bulk fields in the z direction, 



which are holographic images of the hadronic wave func- 
tions. The idea of in-hadron condensates was suggested 



in |22j and developed in [23|, [2J] . Notice that in-hadron 
condensates can be related to vacuum condensates (e.g. 
in the pion case via current algebra) , and therefore there 
is no contradiction or conflict with QCD and chiral per- 
turbation theory 25|. In particular, the main idea of 
Refs. is that the quark and gluon condensates, 

holographically interpreted as spatial effects of bulk fields 
in the z direction, can be interpreted according to the 
AdS/QCD dictionary as finite- volume effects in hadrons. 
As stressed in Ref. [23j the quark and gluon condensates 
have spatial support within hadrons. 

We will show that our approach is consistent with 
model-independent relations and constraints valid in the 
regime of explicitly and spontaneously broken chiral sym- 
metry. The pion is massless only when the chiral sym- 
metry is spontaneously broken (for a vanishing current 
quark mass rh) while it becomes massive via the mech- 
anism of explicit breaking of chiral symmetry encoded 
in the longitudinal part of its LFWF. In addition, we 
generate a finite splitting of the axial- vector and vector 
meson masses and reproduce the Weinberg relations be- 
tween the masses of p(770), ai(I270) and / (600). When 
the current quark masses are included in the formalism 
through the longitudinal LFWF, the pion and other chi- 
ral pseudoscalar fields (kaon and ?y meson) acquire a mass 
according to the Gell-Mann-Okubo-Renner relations. 

We would also like to mention that the classifica- 
tion of the bulk fields propagating in AdSs is similar 
to the classification of hadrons according to the chi- 
ral group [3, l27j . As we stressed before, fields in 
AdSs are classified by unitary, irreducible representa- 
tions of the SO{4,2) group 0,13. The 50(4,2) group 
is decomposed with respect to its maximal subgroup 
SO(4) (g) SU(2), where SO(A) is isomorphic to the group 
SU(2) <g> SU(2). Thus, the bulk fields are character- 
ized by three quantum numbers — minimal energy Eq 
and two spins J\ and J2 — and belong to the rep- 
resentations denoted by D(Eq, J\, J2), which have spe- 
cific chiral properties, because SO(4) is also isomorphic 
to the chiral group SUl{2) ® SUr{2). Because of the 
gauge/gravity duality the energy of the bulk field E is 
identified with A — the dimension of the correspond- 
ing CFT operator. Masses of bulk fields are expressed 
in terms of their energy. In particular, the scalar fields 
S belong to the representation D(Eq,0,Q) with mass 
/i 2 i? 2 = Eq(Eq — 4), which also belongs to the chiral 
representation (0,0) (here R is the AdS radius). There 
are two independent possibilities for a description of vec- 
tor fields: in terms of vectors transforming according 
to the "vectorial" representation D(Eq, i, i) with mass 
/i 2 R 2 = (E — 1)(E — 3), which belong to the chiral 
representation (1,0) © (0,1), and in terms of antisym- 
metric tensors W^ u transforming according to the "ten- 
sorial" representation D(Eq, 1, 0)(BD(Eq, 0, 1) with mass 
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/x 2 i? 2 = [Eq — 2) 2 , which belong to the chiral representa- 
tion (5, |). Finally, the fermions ip with spin 1/2 belong 
to the representation D(E o ,0, ^)®D(E , |,0) with mass 
fiR = E — 2. Fermions belong to the chiral representa- 
tion (0,5) © (5,0). The classification of chiral proper- 
ties of bulk fields in AdS is identical to the classification 



according to a chiral basis [26|, l27|, including two pos- 
sible representations for vector fields, which correspond 
to vector and pseudotensor operators in QCD. In this 
sense, the case of exact chiral symmetry corresponds to 
gauge invariance for massless bulk fields in AdS. The chi- 
ral limit for vector mesons corresponds to conservation of 
the interpolating vector currents on the AdS boundary of 
AdS space. This means that their holographic analogues 
in AdS - vector bulk fields, should be massless or have 
energy Eq = 2 + J = 3 (with J = 1) for both representa- 
tions D(E , ~, 5) axidD(E o ,l,0)®D(E ,0,l). Therefore 
in the case of exact chiral symmetry the twist dimen- 
sion of interpolating operators of vector mesons must be 
A = 2 + J = 3 (with J = 1). However, this is only 
true for a conformal field theory. In order to guarantee 
the correct scaling of hadronic form factors in QCD the 
assignment Eq = A must be connected with the twist 
r of the corresponding interpolating operator expressed 
in terms of L = max \ L Z \ — the maximal value of the z 
component of the quark orbital angular momentum in the 
LF wave function j5i, 28] . In the case of two-parton states 
(mesons) Am = tm = 2 + L, while in the case of three- 
parton states (baryons) A^ =tb + 1/2 = 7/2 + L. This 
means that the scaling of operators in QCD depending 
on L corresponds to the picture of both broken confor- 
mal and chiral symmetries, and hadrons can be classified 
according to the nonrelativistic scheme 2S+1 Lj. Confor- 
mal invariance is broken to the symmetry of the Poincare 
group due to the dilaton field present (living) in the extra 
z dimension, which is dual to the gluon condensate Q- 
The chiral group SU (2) l (g) SU (2)r is broken to the vec- 
tor isospin group SU(2)y due to the coupling of bulk 
fields with the dilaton. 

In a series of papers two versions of soft-wall mod- 
els based on the use of positive 0, H3, [29| and nega- 
tive 13, 3(| dilaton profiles were developed. We showed 
that these models are equivalent in the case of the bound 
state problem due to a specific dilaton field-dependent 
redefinition of the bulk field. It could be shown that 
after such a redefinition the action for the bulk field 
®j = $m 1 ---Mj(%,z) with spin J reads 



S., = / d d xd Zy fg d M t>jd M <t> J - (p 2 , + Vj{z))9j* j ,(1) 



where the AdS metric is specified as ds 2 — e 2A ^ {dx^dx* 1 
A(z) = log(R/z) and R is the AdS 
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radius. Here $j is the symmetric, traceless tensor clas- 
sified by the representation D(Eq, J/2, J/2) with energy 
Eq = A = 2 + L, which is related to the bulk mass 
fij as n 2 R 2 = {E - J)(Eq - 4 + J) = L 2 - (2 - J) 2 . 



This action is most convenient in order to study the 
bound state problem. Versions of the soft-wall model 
with so-called positive or negative dilaton profiles are 
just different manifestations of the action (JT]) after an ap- 
propriate dilaton-dependent redefinition of the bulk field 
— > ^>je ±v( - z ^ 2 . Vj(z) is the effective dilaton poten- 
tial, which has an analytical expression in terms of the 
field ip(z) and the "metric" field A(z) without referring 
to any specific form of their z profiles: 

Vj(z) = e- 2A ^ (y"(z) + (d — 1 — 2 J) <p'(z) A'{z)) . (2) 

The potential Vj(z) breaks both conformal and chiral in- 
variance spontaneously. It was originally suggested in 0] 
to use a quadratic profile for the dilaton field (more pre- 
cisely its z profile) ip(z) = k 2 z 2 , where k is a scale pa- 
rameter, and the conformal metric A{z) = log(R/z) in 
order to obtain Regge behavior for the hadronic mass 
spectra. Note that such a form of the profile arises im- 
mediately if one considers the free action for the dilaton 
in the following form: 



d xdzy/g d M x{x,z)d x(x, z) - H X X (x, z) ,(3) 



where /i 2 = A(A— 4) = —4 is the bulk mass of the dilaton 
with A = 2 (therefore, the dilaton field is the pure scalar 
field). We propose that the Kaluza-Klein expansion for 
the dilaton field is trivial: 



X(x,z) = <p(z) , 



(4) 



where ip(z) is the dilaton z profile, and therefore the dila- 
ton is only living in the z direction. The dilaton should be 
massless (it does not appear in the observable hadronic 
mass spectra). Solving the equation of motion for <p(z) 
with fi 2 , R 2 = —4 and A(z) = \og[R/z) and assuming the 
power behavior of y^jWe find (p(z) ~ z 2 which confirms 
the conjecture of Ref. [4|. 

As we stressed before, the quadratic dilaton can be 
considered as the expectation value of the scalar bulk 
field with dimension 2, which is holographically dual to 
the dimension-2 gluon operator (^4 2 jlnin ) [llj — a nonlo- 
cal gauge-invariant gluon condensate coinciding with the 
gauge-noninvariant gluon condensate of the same dimen- 
sion (A 2 ) in the Landau gauge. The dimension-2 gluon 
condensate has been discussed in the literature (see e.g. 
Refs. [1,[21, 31 1). Note, the interpretation of the dilaton 
as the quantity dual to the condensate of the dimension-2 
operator has been given in the framework of the soft-wall 
model where the dilaton was introduced in the warp- 
ing factor, breaking the conformal-invariant background 
metric. Inclusion of a more complicated form of the dila- 
ton potential (e.g. taking into account self-interaction 
terms) can be viewed as a further extension of the soft- 
wall models based on the quadratic dilaton profile. 

Notice that a quadratic form of the z profile of the 
dilaton is not unique. For example, in the Liu-Tseytlin 
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model (a type of top-down AdS/QCD approach) [32j the 
conformal invariance was violated by the dilaton taken in 
the form = 1 + qz 4 , where the parameter q, accord- 
ing to the AdS/QCD dictionary [33(, was related to the 
matrix element of a QCD operator: the scalar (tv{G 2 )) 
and pseudoscalar (tr(G M „G Mt/ )) gluon condensates. On 
the other hand, the scalar dimension-4 gluon condensate 
is connected to the quark vacuum condensate via the de- 
coupling relation [34j 

1 



(Q|gg|0) 



12m 



R tr(G^)|0) 



(5) 



derived in the leading order of a 1/m expansion, where 
m is the mass of a heavy quark or the constituent 
mass in the case of light quarks. Therefore, nonzero 
dimension-4 gluon condensates signal the existence of 
nonzero dimension-6 quark condensates, which are man- 
ifestations of the spontaneous breaking of chiral symme- 
try. Therefore we see that the breaking of conformal 
symmetry is connected with the breaking of chiral sym- 
metry. Actually, this idea is not new (see e.g. discussion 
in Ref. (H). 

Next we explain why the dilaton is also responsible for 
the spontaneous breaking of chiral symmetry. With the 
quadratic z profile of the dilaton field <p(z) = k 2 z 2 the 



meson spectrum is given by the master formula 13 1 

L + J 



4k< 



(6) 



where k is the dilaton parameter, which signals the spon- 
taneous breaking of conformal and chiral symmetry and 
leads to a discrete mass spectrum of hadrons. An addi- 
tional mechanism of spontaneous breaking of chiral sym- 
metry is encoded in the bulk mass p: 2 , which explicitly 
depends on L and forbids parity doubling (e.g. between 
vector and axial mesons). The corresponding breaking 
term in the mass formula is 5M 2 — 2k 2 (L— J). It means 
that we can interpret Eq. ([6]) as 



M 



M 2 + 8M 2 



(7) 



where M 2 = An 2 (n + J) is the term corresponding to the 
parity doubling limit. 

The hadronic wave functions are identified with the 
profiles of AdS modes $„(z) in the z direction: 



= V r(n + L + l) KL+lz2+Le ' K2Z2/2£ " (K2z2),(8) 

which have the correct behavior in both the ultraviolet 
and infrared limits, 

$„l(z) z 2+L at small z, <& n L{z) -> at large z (9) 

and are normalized according to the condition 

oo oo 

-/: . ., / , . (1Q) 



-$Kl{z) = / dz<t>i L {z) = 1, 



where $„t(z) = z 3 / 2 0„l(z). At z — > the scaling of 
the bulk profile is identified with the scaling of the cor- 
responding mesonic interpolating operator t = 2 + L. 
As we mentioned in the Introduction, r depends on L 
(instead of J as in CFT) because we model QCD and 
should reproduce the scaling of hadronic form factors. As 
we stressed before, dependence on L means spontaneous 
breaking of chiral invariance, which is expected because 
after the introduction of the dilaton field we broke the 
conformal or gauge invariance acting in AdS space. As 
we noted before, the chiral group is isomorphic to the 
subgroup of S"0(4, 2). 

Let us discuss the consistency of the mass formula ^ . 
First of all, we reproduce the massless pion M 2 = for 
n = L = J = 0, which is consistent with the picture of 
spontaneous breaking of chiral symmetry. Also, at z = 
all bulk profiles vanish, which corresponds to pointlikc 
hadrons (zero scale limit) and we restore conformal, chi- 
ral, and gauge invariance associated with the AdS group 
50(4, 2). In this sense we do not have any contradictions 
with chiral invariance as mentioned in Ref. [27j |. because 
for z — it is restored for harmonics with any L and for 
finite z we live in the phase of spontaneously broken chi- 
ral invariance (see further details in (36|). For finite z the 
profile for the ground state p meson is defined by their 
leading twist corresponding to L — 0, which is consistent 
with the statement of Ref. [27j , that the ground state p 
meson is almost a pure S wave in the infrared. On the 
other hand, the bulk profile dual to the ground state p 
meson is defined in the extra dimension. Therefore it has 
no correspondence to the chiral-invariant superposition 
of Fock states with L = (S'-wave) and L = 2 (D-wave) 
discussed in (2^]. Note that both the electromagnetic 
field and J = L = 1 vector mesons have the same inter- 
polating twist-3 operator. They also have the same dual 
field in AdS space — a massless vector bulk field and 
the bulk-to-boundary propagator of the electromagnetic 
field in the timelike region [28|, |3tJ , which has the poles 
corresponding to the excited states with J = L = 1: 



P 



M 



4k 2 (ti 



1 



(11) 



as expected in the chiral-invariant limit. 

Let us mention other interesting results deduced from 
the meson mass formula ([6]). There are two relations 
between the masses of the ground states of vector p(770), 
axial-vector ai(1270), and scalar /o(600) mesons, 



M V2 = 2k : 



M 



fo 



V2k (12) 



which are consistent with the predictions done in Ref. 
in the same limit when chiral symmetry is spontaneously 
broken. Also, the vector and axial- vector multiplets are 
not degenerate in our approach (even for higher values 
of J), because of the finite mass splitting of axial- vector 
and vector mesons states: 



Mi 



M 2 



2k 2 



(13) 
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This means that we do not have parity doubling. An 
amazing fact is that the same prediction relating the 
masses of p and a x mesons was also obtained [391 ] on 
the basis of spectral function sum rules at M n = and 
using the Kawarabayashi-Suzuki-Riazuddin-Fayyazuddin 
formula (4pj for the pirn coupling. 

For the value k — 500 MeV used in this paper the 
masses of the p and a\ mesons are well reproduced in 
comparison to data 

M p = 721 MeV, (data : 775.49 ± 0.34 MeV) , 
M ai = 1010 MeV , (data : 1230 ± 40 MeV) , (14) 

while our result for the /o mass should be considered 
as a prediction. We get Mf — 721 MeV, which per- 
fectly agrees with a model-independent result based on 

M f0 = 



analyticity and unitarity of the S matrix 41 



735.0 ±6.1 MeV. But also note that a recent compila- 
tion [42| indicates a mass of Mf — 446 ± 6 MeV as 
deduced from the pole position in the process amplitude. 

Next we consider the inclusion of the longitudinal part 
of the LFWF. We will show that this extension is im- 
portant to include in our formalism the dependence on 
the current quark masses and therefore explicit break- 
ing of chiral symmetry. As we stressed before, the first 
step in this direction was done in Refs. 10]. The au- 
thors proposed to write down the mesonic two-parton 
wave function in a factorized form, as a product of trans- 
verse </vil(C): longitudinal f(x, mi, 777,2) and angular e tm ^ 
modes. In Ref. [13( we further proposed to do such sepa- 
ration in a more convenient form, factorizing in addition 
\J x(\ — x) — the Jacobian of the £ — > |bj_| coordinate 
transformation: 



f(x,mi,m 2 ) 




(15) 



In Ref. 13] we used a Gaussian ansatz for the longitu- 
dinal part of the LFWF and treated the quark masses 
mi and 777,2 as constituent quark masses. Here we fol- 
low Refs. [lH [r| and consider current quark masses, 
and by an appropriate choice of the longitudinal wave 
function f(x, mi, 7712) we get consistency with QCD in 
both sectors of light and heavy quarks. We generate the 
masses of light pseudoscalar mesons in agreement with 
the picture resulting from explicit breaking of chiral sym- 
metry — in the leading order of the chiral expansion 
the masses of pseudoscalar mesons are linear in the cur- 
rent quark mass. In this vein we also guarantee that the 
pseudoscalar meson masses satisfy the Gell-Mann-Oakes- 
Renner (GMOR) relation for the pion mass 



Ml = 2m B 



(16) 



and the Gell-Mann-Okubo (GMO) relation between the 
masses of pion, kaon, and r\ meson: 



4M% 



Mi 



3Mi 



(17) 



In previous equations m = (m u + md)/2 is the average 
mass of u and d quarks (we work in the isospin limit 
m u = rn d ), B = |(0|uu|0)|/F^ is the quark condensate 
parameter, and F„ is the leptonic decay constant. Note, 
the condensate parameter B is related to the coupling 
constant of the pseudoscalar density to the pion G^ [25 1 
(or in-hadron condensate of pion 2^, 24 1): 



B 



2K 



(0\qi-f 5 T i q\TT k 



6 ik Gir 



(18) 



In the sector of heavy quarks we get agreement with 
heavy quark effective theory and potential models for 
heavy quarkonia. In the heavy quark mass limit vtlq — ¥ 
00 we obtain the correct scaling of the leptonic decay 
constants for both heavy- light mesons /q<j ~ 1/ ^Jvtiq 
and heavy quarkonia /qq ~ s/tuq and f c i ~ m c / ^rrib 
at m c -C mf,. In this limit we also generate the correct 
expansion of heavy meson masses 



M Qq = 
Man = 



m Q H 
2m Q 



A 



E 



0(l/mq), 
+ 0(l/m Q ), 



(19) 



where A is the approximate difference between the masses 
of the heavy-light meson and the heavy quark, E is the 
binding energy in heavy quarkonia, and their splittings, 
e.g. between vector and pseudoscalar states of heavy- 
light mesons, become 



1 

m Q 



(20) 



We choose the longitudinal wave function in the form 

f(x,m u m 2 ) = Nx ai {l~x) a2 (21) 

where N is the normalization constant fixed from 
1 



1 = J dx f (x, mi, 7712) 
" 



(22) 



and ai, Oi2 are parameters that will be fixed in order to 
get consistency with QCD. 

In the present paper the physical quantities of interest 
are the mass spectrum M^j and lepton decay constants 
fni of mesons, which are given by the expressions fl3| 



L + J 



x 1 — x 



f (x,mi,m 2 ) , (23) 



V6 



f M = k — / dxy/x(l -x)f(x, mi,m 2 ). (24) 
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Using our ansatz for the longitudinal wave function ([21] 
we calculate the leptonic decay constant and the correc 
tion to the mass spectrum analytically 13j: 

Ml, = W L+J 



(l + 2ai 
V6 r( 



2a 2 ) ^ + ^ 



fia = k — 



2a 



| + « 1 )r(f 



2a; 



■ a 2 J 



T(3 + ai + a 2 ) 

/ T(2 + 2ai + 2a 2 ) 
r(l + 2ai)r(l + 2a 2 ) 



(25) 



Next we analyze different types of mesons. We start with 
the light pseudoscalar mesons. Here we want to incorpo- 
rate a mechanism for explicit breaking of chiral symmetry 
in order to reproduce the GMOR and GMO relations in 
the leading order of the chiral expansion. For this pur- 
pose we set ai — mj/(2B), which is consistent with the 
ideas and results of Refs. [TBI, [l6| and which means that 
ai vanishes in the chiral limit m t — > 0. The leptonic de- 
cay constants of both the pion and kaon are finite and 
degenerate in the chiral limit 0, UH : 

f P = ^- (26) 

Concerning the longitudinal wave functions of other light 
mesons, the choice of the a; parameters does not nec- 
essarily coincide with the one used for the pseudoscalar 
mesons. But for simplicity we will use the same universal 
longitudinal wave function for all light mesons. In this 
case the leptonic decay constants of it, K, and p° mesons 
are degenerate and the corresponding constants for u and 
<fi mesons are related via the SU(3) flavor conditions: 

,■ / _ 3 /0 ^6 
Jp 



fir I j 



K 



3/. = ^ = k ^ ■ (27) 

Next we consider the heavy-light mesons. Here we should 
reproduce the heavy quark mass expansion of the heavy- 
light meson masses, the mass splitting of vector and pseu- 
doscalar states, and the l/y'mg scaling of the leptonic 
decay constants. All these constraints are fulfilled if the 
ai parameters are fixed as aQ — a, & flavor-independent 
constant, while the parameter a q 

a q = — ^ 1 
m Q V 



must be 
2m Q J 2 



(28) 



For this choice the results for the mass spectrum and 
leptonic decay constants, in the leading order of the chiral 
expansion, and the leading and next-to-leading order of 
the heavy quark mass expansion, are 

L + ,r 



fQq 




(29) 



— ~\— ■ (30) 



Finally, for heavy quarkonia we fix the parameters aQ i 
as 



a Qz = -rh- I 1 



E 



+ 0(—) (31) 



AE \ 2(m Ql +TOQ 2 )/ \m Qi 
and then we get the following results for the spectrum 
L + J s 



Q1Q2 



+ (m Ql + m Q2 + Ef . (32) 



It was shown in [43J that the trajectories of bottomia 
states deviate from linearity as we also discussed in fl3j j . 
This effect can be related to the one-gluon exchange 
term, which results in an additional Coulomb-like inter- 
action between quarks V(r) = — 4a s /3r, where a s is the 
strong coupling constant. Its contribution to the mass 
spectrum M is negative and proportional to the quark 
mass squared [43| . For light and heavy- light mesons this 
term can be safely neglected, while this is not the case 
for heavy quarkonia (especially for bottomia states). Ex- 
tending the results of Refs. |43| to the general case of a 
meson containing constituent quarks with masses m,Q 1 
and niQ 2 , we get the following expression for the shift of 
M 2 due to the color Coulomb potential: 



am; 



64a 2 s m Ql mQ 2 
~9(n + L + If 



(33) 



where a s is the strong coupling considered as a free pa- 
rameter. Therefore, the final expression for the heavy 
quarkonia spectra is given by the master formula: 



Mi 



An 2 [ n + ) + (m Ql + m Q . 2 + E) 2 



64a 2 m Ql m Q2 
%{n + L + \f 



(34) 



For the leptonic decay constants of unflavored quarkonia 
we get the following result in leading order of the 1 /uiq 
expansion: 



(2^)3/4 [ mQ J 
and for the decay constant of B c meson: 
2kV6 (m c /E)V 4 



fcb — 



TT 3 / 4 (m b /E) ' 



(35) 



(36) 



In the case of the B c meson we additionally apply the 
condition m c <C mj. Equations (|3"5|) and (|36p can be 
combined into a general formula for the leptonic decay of 
heavy quarkonia: 



/Qi< 



2^76 (fi QlQ2 /Ef^ 
TT 3 / 4 (M Q /E) 



(37) 



where ^q 1 q 2 — mQ 1 mQ 2 /(m,Q 1 + mQ 2 ) is the reduced 
mass of heavy quarkonia and Mq = vtiq 1 + mQ 2 . In 
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order to get the correct scaling [44[ of the leptonic decay 
constants of heavy quarkonia we propose the following: 
in the case of heavy quarkonia the dilaton parameter n 
should be flavor dependent and scale as Pq^q 2 ' M^ 2 . 
Note that different dilaton parameters for light mesons 
and heavy quarkonia were also considered before (see e.g. 
Ref. 45]). Here we use the following ansatz for k: 



K = /3 



MQ1Q2 
E 



1/4 



Mo 
E 



1/2 



(38) 



where /3 = O(mg). For unflavored quarkonia the final 
result for the leptonic decay constant reads: 



_ V3 lm c , 
JQQ — P ~T7I \ "FT ■ 



For the B c meson we get 

V6 m c /E 



fob = 2/3 



7T 3 / 4 ^/E 



(39) 



(40) 



Finally we present the numerical results. The parameters 
are fixed to the following values. For the current quark 
masses we use 

m u = md = rh = 7 McV , 
m a = 24m = 168 MeV , 

m c = 1.275 GeV , m b = 4.18 GeV . (41) 

For the strong coupling constants a s we use the following 
set of parameters: 

a s (cc) = 0.45, a s (cb) = 0.383, a s (&&) = 0.27. (42) 

Note that for the light and heavy-light mesons we use the 
universal parameter k = 500 MeV. The parameters A for 
heavy-light mesons are fixed from the mass difference of 
the experimental values of ground state of D, D s , B, B s 
mesons and corresponding heavy quark mass as 

A cq = 0.595 GeV , A cs = 0.695 GeV , 

A bq = 1.1 GeV , A &s = 1.19 GeV . (43) 

The binding energies for heavy quarkonia are fixed as 

E cc = 0.795 GeV , E cb = 1.25 GeV , 



E bb = 1.45 GeV. 



(44) 



The set of /3 couplings, defining the heavy flavor depen- 
dence of the dilaton parameter k, is fixed in the case of 
heavy quarkonia as 

/3(cc) = 0.36 GeV , (3(cb) = 0.32 GeV , 

/3(bb) = 0.41 GeV . (45) 



(3 = 0(itiq). In Tables Hl lVIl we present the numerical 
results both for the mass spectrum and the leptonic de- 
cay constants of light, heavy- light and heavy quarkonia 
in leading order of the chiral and heavy-quark mass ex- 
pansion. Results are for most of the part in reasonable 
agreement with data. We would like to stress that we sig- 
nificantly improved the description of mesonic properties 
in comparison to our previous efforts (l3 |. 

In conclusion, we demonstrated that in the soft-wall 
model, where conformal or S*0(4, 2) gauge invariance is 
broken the same applies for chiral invariance which is also 
broken (spontaneously). This is manifested in the L de- 
pendence of twists of interpolating operators of hadrons 
and in their observables such as mass spectra and form 
factors. In the limit of chiral invariance one could ex- 
pect a dependence on J. However, the J dependence 
is specific for conformal field theories and not for QCD. 
Namely, the L dependence is dictated by the scaling of 
hadronic form factors at higher Q 2 .A restriction to a spe- 
cific L (in our case to the minimal L for a specific hadron) 
is a manifestation of spontaneous breaking of chiral sym- 
metry. Chiral symmetry is restored in the exact limit 
z = (not for small z), which is trivial and corresponds to 
the restoration of chiral symmetry for all hadrons having 
components with an adjustable value of L. We demon- 
strated explicitly that the present approach is consistent 
with model-independent predictions obtained in the case 
of spontaneous broken chiral symmetry: i) a massless 
pion, ii) the Weinberg sum rule relating masses of p, a\ 
and /o mesons. In the chiral limit the hadron eigenstates 
are superpositions of components with different values of 
L. Finally, we demonstrated how to consistently con- 
struct the longitudinal LFWF and include the current 
quark mass dependence. The latter is consistent, in the 
light quark sector, with the mechanism of explicit break- 
ing of chiral symmetry, and in the heavy quark sector 
with the heavy quark spin-flavor symmetry. 
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These parameters are nearly the same for all quarko- 
nia states, which is consistent with proposed scaling 
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TABLE I: Masses of light mesons. 



Meson 


n 


r 

L 


o 
o 


Mass 


MeV] 


TV 


0,1,2,3 


U 


U 


Wtt(140) = 140 


71 /T 1 A 1 A 

M w( i300) = 1010 


71 /T 1 A 1 

A^(i800) = 1421 


71 /T 1 TO O 

A«7r(4 s ) = 1738 


K 


U 


A 1 O O 


U 


Mk = 495 


Mffi(1270) = 1H6 


71 /f 1 /I AO 

Mk 2 (i77o) = 1498 


71 X 1 OA 1 

Mk 3 = 1801 


V 


A 1 O Q 

U,l,z,o 


U 


U 


M n(ls) = 566 


71 /T 1 1 /I A 

M^ (2s ) = 1149 


71 /T i roo 

M v ( 3a ) = 1523 


71 /T 1 OOO 

M v ( As ) = 1822 


Jo [nn\ 


A 1 1 O 


1 


1 


71 /f 701 

m /o(1p) - 721 


7l/f 1 OQQ 

M /o(2 p) - 1233 


7l/f 1 

m /o(3p) - 1587 


m /o(4p) - 1S76 


Jo [ss\ 


n i o q 
U,l,z,o 


1 


1 


M Mlp) = 985 


M M2p) = 1404 


M M 3 P ) = 1723 


M M4p) = 1993 


p(770) 


0,1,2,3 





i 


M p(770) = 721 


A/p(i450) = 1233 


M p{1700) = 1587 


M p(4s) = 1876 


cj(782) 


0,1,2,3 





i 


M w(782) = 721 


^w(1420) — 1233 


M u(l650) = 1587 


M w(4s) = 1876 


0(1020) 


0,1,2,3 





i 


M^ (ls) = 985 


M H2a) = 1404 


M^ (3s) = 1723 


M H4s) = 1993 


ai(1260) 


0,1,2,3 


1 


i 


M ol( i p) = 1010 


M ol(2p) = 1421 


M ol(3p ) = 1738 


M ol(4p ) = 2005 



TABLE II: Masses of heavy-light mesons. 



Meson 


J p 


n 


L 


S 


Mass 


MeV 




D(1870) 


0" 





0,1,2,3 





1870 


2000 


2121 


2235 


L>*(2010) 


1" 





0,1,2,3 


1 


2000 


2121 


2235 


2345 


D s (1969) 


0" 





0,1,2,3 





1970 


2093 


2209 


2320 


D s *(2107) 


1" 





0,1,2,3 


1 


2093 


2209 


2320 


2425 


B(5279) 


0" 





0,1,2,3 





5280 


5327 


5374 


5420 


B*(5325) 


1" 





0,1,2,3 


1 


5336 


5374 


5420 


5466 


B s (5366) 


0" 





0,1,2,3 





5370 


5416 


5462 


5508 


B 3 *(5413) 


1" 





0,1,2,3 


1 


5416 


5462 


5508 


5553 



TABLE III: Masses of heavy quarkonia cc, bb and cb 



Meson 


J p 


n 


L 


s 


Mass 


MeV] 


r/ c (2980) 


0" 


0,1,2,3 








2975 


3477 


3729 


3938 


V>(3097) 


1" 


0,1,2,3 





1 


3097 


3583 


3828 


4032 


Xco(3415) 


0+ 


0,1,2,3 


1 


1 


3369 


3628 


3843 


4038 


Xd(3510) 


1+ 


0,1,2,3 


1 


1 


3477 


3729 


3938 


4129 


Xc2(3555) 


2+ 


0,1,2,3 


1 


1 


3583 


3828 


4032 


4219 


r/ 6 (9390) 


0" 


0,1,2,3 








9337 


9931 


10224 


10471 


T(9460) 


1" 


0,1,2,3 





1 


9460 


10048 


10338 


10581 


XM)(9860) 


0+ 


0,1,2,3 


1 


1 


9813 


10110 


10359 


10591 


X6i(9893) 


1+ 


0,1,2,3 


1 


1 


9931 


10224 


10471 


10700 


X62(9912) 


2+ 


0,1,2,3 


1 


1 


10048 


10338 


10581 


10808 


B c (6277) 


0" 


0,1,2,3 








6277 


6719 


6892 


7025 



TABLE IV: Decay constants fp of pseudoscalar mesons in MeV. 





Data [42l 


Our 


TV 


130.4 ± 0.03 ± 0.2 


153 


K~ 


156.1 ± 0.2 ± 0.8 


153 


D+ 


206.7 ±8.9 


207 


Dt 


257.5 ±6.1 


224 


B~ 


193 ± 11 


163 


B° 


253 ± 8 ± 7 


170 


B c 


489 ± 5 ± 3 [46] 


489 



TABLE V: Decay constants fy of vector mesons with open flavor in MeV. 



Meson 


Data 


Our 


P + 


210.5 ±0.6 [42] 


216 


D* 


245 ± 20i^ [4JJ 


207 


Dt 


272 ± 16±1 [48] 


224 


B* 


196 ± 2A±f [4JJ 


163 


Bt 


229 ± 20t*g [47] 


170 



TABLE VI: Decay constants /y of vector mesons with hidden flavor in MeV. 



Meson 


Data [42] 


Our 


P° 


154.7 ± 0.7 


153 


UJ 


45.8 ± 0.8 


51 


4> 


76 ± 1.2 


72 


J/4, 


277.6 ± 4 


223 


T(ls) 


238.5 ± 5.5 


170 



